Discrete breathers in nonlinear lattices: Experimental detection in a Josephson array 
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We present an experimental study of discrete breathers in an underdamped Josephson-j unction 
array. Breathers exist under a range of dc current biases and temperatures, and are detected by 
measuring dc voltages. We find the maximum allowable bias current for the breather is proportional 
to the array depinning current while the minimum current seems to be related to a junction retrap- 
ping mechanism. We have observed that this latter instability leads to the formation of multi-site 
breather states in the array. We have also studied the domain of existence of the breather at different 
values of the array parameters by varying the temperature. 



Discrete breathers are a new type of excitation in non- 
linear lattices. They are characterized by an exponential 
localization of the energy. This localization does not oc- 
cur in linear systems and it is different from Anderson 
localization, which is due to the presence of impurities. 
Thus, discrete breathers are also known as intrinsic lo- 
calized modes. 

Breathers have been proven to be generic solutions 
for the dynamics of nonlinear coupled oscillator systems 
jl],D by the use of the novel mathematical technique of 
the anti-integrable limit ||. They have been extensively 
studied ^-||] and have been proposed to theoretically 
exist in diverse systems such as in spin wave modes of 
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antiferromagnets |9|], DNA denaturation [10], and the 
dynamics of Josephson-junction networks Also, 
they have been shown to be important in the dynam- 
ics of mechanical engineering systems |l3|,ft4j ■ Although 
a number of experiments have been proposed, discrete 
breathers have yet to be experimentally generated and 
measured. 

In this Letter, we present, to our knowledge, the first 
experimental study of discrete breathers in a spatially ex- 
tended system. We have designed and fabricated an un- 
derdamped Josephson-junction ladder which allows for 
the existence of breathers when biased by dc external 
currents. We have developed a method for exciting 
breathers and explored their existence domain and in- 
stability mechanisms with respects to the junction pa- 
rameters and the applied current. 

A Josephson junction consists of two superconduct- 
ing leads separated by a thin insulating barrier. Due to 
the Josephson effect, it behaves as a solid-state nonlinear 
oscillator and is usually modeled by the same dynami- 
cal equations that govern the motion of a driven pendu- 
lum p5[ , p!6[ : i = ip + Tip + sin ip. The response of the 
junction to a current is measured by the voltage of the 
junction which is given by v = {^a/2-K)dip / dt. By cou- 
pling junctions it is possible to construct solid-state phys- 
ical realizations of different models such as the Frenkel- 
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FIG. 1. Current-voltage characteristic of an anisotropic 
Josephson junction array. The hysteresis between the de- 
pinning current (Ide P ~ 2 mA) and the retrapping current 
(I r ~ 0.2 mA) is shown. Inset: Schematic of the anisotropic 
ladder array. Vertical junctions have four times the area of the 
horizontal ones. Current is injected through bias resistors, Rb, 
of 25 Q. These resistors are larger than the vertical junction's 
normal state resistance of 5 Q to assure a more homogeneous 
current drive and are also larger then the horizontal's normal 
state resistance of 20 Q to minimize effects on their dynam- 
ics. There are voltage probes in the fourth, fifth, sixth and 
ninth vertical junctions to measure Va, V5, Ve and Vg. The 
voltage probes can also be used to measure the top horizontal 
junctions in the middle which we denote as Vat and Vst- 



Kontorova p7| model for nonlinear dynamics and the 2D 
XY model ]18| ] for phase transitions in condensed matter. 
Moreover, since the parameters, such as T(T), vary with 
temperature, a range of parameter space can be studied 
easily with each sample. 

The inset of Fig. [j] shows a schematic of the anisotropic 
ladder array. The junctions are fabricated using a Nb- 
AkjOz-Nb tri-layer technology with a critical current den- 
sity of 1000 A/cm 2 . The current is injected and extracted 
through bias resistors in order to distribute the current 
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as uniformly as possible through the array. These re- 
sistors are large enough so as to minimize any delete- 
rious effects on the dynamics. The anisotropy of the 
array is defined by h as the ratio of areas of the hori- 
zontal to vertical junctions. In our arrays h — 1/4 and 
h = I ch /I cv = R v /R h = C h /C Vl and r„ = T h = T. As 
shown in the schematic, we have placed voltage probes 
at various junctions in order to measure the voltages of 
both horizontal and vertical junctions. 

In Fig. [l] we show a typical current- voltage, IV, charac- 
teristic of the array. As the applied current increases from 
zero we measure the average voltage of the 9-th junction. 
The junction starts at a zero-voltage state and remains 
there until it reaches the array's depinning current Idep 
at about 2 mA. The depinning current can roughly be 
understood as the sum of the vertical junction's intrin- 
sic critical current I cv and the small circulating Meissncr 
current around the array. In a pendulum analogy, the 
critical current is equivalent to the critical torque, which 
is just sufficiently strong enough to force the pendulum 
to start rotating. When the current is larger the junction 
switches from zero-voltage state to the junction's super- 
conducting gap voltage, V g , which at this temperature 
is 2.5 mV. At this point all of the vertical junctions are 
said to be rotating and the array is in its "whirling state" . 
One of the effects of this gap voltage is to substantially 
affect the junction's resistance, and thereby damping, in 
a complicated nonlinear way. The current can be further 
increased until the junction reaches its normal state and 
it behaves as a resistor, R n , of 5 O. As the current de- 
creases the junction returns to the gap voltage and then 
to its zero- voltage state at the retrapping current, I r , 
of 0.2 mA. The hysteresis loop between Id ep and is 
due to our underdamped junctions: the inertia causes the 
junctions to continue to rotate when the applied current 
is lowered from above its critical value. 

It is this hysteresis loop that allows for the existence 
of breathers in the ladder with dc bias current. In this 
current range the zero-voltage (V = 0) and rotating 
(V = V g ) solutions coexist. Then, a discrete breather 
in the ladder corresponds to when one vertical junction 
is rotating while the other vertical junctions librate. This 
solution is easy to conceive in the limit where the vertical 
junctions are imagined to be completely decoupled. How- 
ever, whether a localized solution can exist in the ladder 
will be determined by the strength of the spatial coupling 
between vertical junctions. This coupling occurs through 
three mechanisms: flux quantization, self and mutual in- 
ductances of the meshes, and the horizontal junctions. 
Though the effective coupling is a complicated function 
of the array parameters, it is most strongly controlled 
by h. If the anisotropy h is too large, then the vertical 
junctions will not support localized solutions that can be 
excited by dc currents. It has been determined from sim- 
ulations of the system jl2| that h — 1/4 will allow for the 
existence of breathers in our ladders. 
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FIG. 2. Schematic of solutions in the array: (a) the 
whirling state where vertical junctions are rotating and hori- 
zontal ones librate; (b) the zero-voltage state where there are 
no rotating junctions; (c) a breather solution where the fifth 
vertical junction and the nearest horizontal neighbors rotate 
and the other vertical and horizontal junctions librate; (d) a 
multi-site breather solution where two vertical junctions ro- 
tate. 
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FIG. 3. Measured averaged voltages of five junctions in the 
center of array as the applied current is increased. We have bi- 
ased the ladder at 1.4 mA and excited a breather as indicated 
in the text. Then the applied current is increased. There are 
two regions in the IV. Below ~ 2 mA we see the breather; 
above, the breather becomes unstable and the array switches 
to the whirling state. 

Figure || shows some possible solutions for the states 
of our ladder. Graph (a) is the whirling state with ev- 
ery vertical junction rotating as indicated by the arrows. 
This is the state when all the junctions of the array have 
switched to the gap voltage. Graph (b) shows the zero 
voltage state with no rotating junctions. Graph (c) de- 
picts a single-site breather solution. Here, one of the 
vertical and the horizontal neighboring junctions rotate. 
The horizontal junctions allow the vertical junction to 
rotate with a mean voltage without an overall increase 
of the stored magnetic energy. Graph (d) shows a two- 
sited breather where two vertical junctions rotate. We 
have experimentally detected these types of localized so- 
lutions (c and d) by measuring the average dc voltage of 
the junctions as labeled in (c). 

For our experiments, we have developed a simple re- 
producible method of exciting a breather: (i) bias the 
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array uniformly to a current below depinning current; 
(ii) increase the current injected into the middle vertical 
junction [labeled V5 in Fig. 0(c)] until its voltage switches 
to the gap; (iii) reduce this extra current in the middle 
junction to zero. 

For example, to prepare the initial state in Fig. ^ we 
started by increasing the applied current to 1.4 mA which 
is below Idep- At this point the array is in the zero- voltage 
state. We then add an extra bias current to the middle 
junction (number 5) until it switches to the gap voltage 
of 2.5 mV and then we reduce this extra bias to zero. In 
a sense we have prepared the initial conditions for the 
experiment. We can now increase the uniform applied 
current while simultaneously measuring the voltages of 
the vertical junctions (V4, V5 and Vq) and the top two 
horizontal junctions, Vat and V5T, as labeled in Fig. ||(c). 

Figure || shows the result after we have excited the 
breather and we have increased the array current. Close 
to the initial current of 1.4 mA only the fifth vertical 
junction is at V g and both the fourth and sixth verti- 
cal junctions are in the zero-voltage state. This is the 
breather state shown in Fig. ||(c) and in essence the signa- 
ture of the localized breather: a vertical junction is rotat- 
ing while its neighboring vertical junctions do not rotate. 
We also see that both neighboring horizontal junctions 
have a voltage magnitude that is precisely half of this 
value (V iT = -V 5 t = V 5 /2 and V 4 = V 6 = 0). Both the 
magnitude and the sign can be understood by applying 
Kirchoff's voltage law to top-bottom voltage-symmetric 
solutions, as sketched in Fig. |^(c). The voltage of the 
top horizontal junction is equal to the negative of the 
bottom one. Since the voltage drops around the loop 
must be zero, the horizontal voltages must be half that 
of the active vertical junction voltage. As we increase 
the current the breather continues to exist until the ap- 
plied current approaches i+ « 2 mA. At this point the 
horizontal junctions switch to a zero-voltage state while 
all of the vertical junctions switch to V g — 2.5 mV. The 
array is now in its whirling state as drawn in Fig. |^(a) 
where V4 = V e — V 5 = V g and V4T = V 5 t = 0. 

If we excite the breather again but instead of increas- 
ing the applied current we decrease it, we measure curves 
typical of Fig. ||. As explained above, we prepare the 
array in an initial condition with a breather located in 
junction 5 at 1.4 mA. We then decrease the applied cur- 
rent slowly. We start with the signature measurement 
of the breather: junction five is rotating at V g while V4 
and Ve = 0. We also see that the horizontal junctions 
have the expected value of V g /2. As the current is de- 
creased the breather persists until the array is biased at 
0.8 mA. The fourth vertical junction then switches to the 
gap voltage while V±t switches to a zero voltage state. 
The resulting array state is sketched in Fig. ||(d) with 
Vi = V 5 = V g while V 5T = Vg/2 and V 4T = V 6 = 0. 
The single-site breather has destabilized by creating a 
two-site breather. 
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FIG. 4. Measured average voltages of five junctions in the 
center of array as the applied current is decreased. We have 
biased the ladder at 1.4 mA and excited a breather as in- 
dicated in the text. Then the applied current is decreased. 
There are three regions in the IV. The breather state be- 
comes unstable at around 0.8 mA leading to a sequence of 
shifts at the gap voltage that are interpreted as multi-site 
breather states. The array reaches its zero-voltage state at 
about 0.2 mA. 

As the applied current is further decreased beyond the 
single-site breather instability at « 0.8 mA, the voltage of 
the fourth and fifth vertical junctions decreases but then 
suddenly jumps back to . Then the voltage decreases 
again, and it again jumps back to V g . This second shift 
corresponds to the sixth junction switching from the zero 
voltage state to the gap voltage. At this current bias, 
all of the three measured vertical voltages are rotating. 
There is a further jump of the voltage as the current de- 
creases. Finally, at 0.2 mA all of the vertical junctions 
return to their zero- voltage state via a retrapping mech- 
anism analogous to that of a single pendulum. 

From these experiments and corroborating numerical 
simulations we speculate that this shifting of the volt- 
age back to V g corresponds to at least one vertical junc- 
tion switching from the zero-voltage state to the rotat- 
ing state. The shapes of the IV curves in this multi-site 
breather regime are influenced by the redistribution of 
current when each vertical junction switches. This re- 
distribution may also govern the evolution of the system 
after each transition to one of the other possible breather 
attractors in the phase space of the array. However, the 
exact nature of the selection process is not yet under- 
stood. 

The above data was taken at a temperature of 5.2 K. 
We found four current values of importance: the cur- 
rent when the array returns to the zero- voltage, I r ; the 
maximum zero- voltage state current, Idep] the maximum 
current the breather supports, i+; and the minimum cur- 
rent I- . By sweeping the temperature we can study how 
the current range in which our breather exists is affected 
by a change of the array parameters. 
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FIG. 5. Existence region of the breather in the current-P 
plane. T was varied by changing temperature. 1+ is the max- 
imum applied current the breather supports, while I- is the 
minimum current. 

Figure || shows the results of plotting the four special 
current values versus T. To calculate how the junction 
parameters vary with temperature we take I cv (0)R n = 
1.9 mV and assume that the junction critical current fol- 
lows the standard Ambegaokar-Baratoff dependence [^9) . 
We estimate T from I r by the relation I r /NI cv = (4/7r)r 
, where N is the number of vertical junctions in the 
array. The other relevant parameter is the dimensionless 
penetration depth, Aj^ = <$>o/2irL s I cv , which measures 
the inductive coupling in the array. The loop induc- 
tance L s is estimated from numerical modeling of the 
circuit. By changing the temperature of the sample, we 
vary the I cv of the junction and hence change T and Aj_. 
In this sample, the junction parameters can range from 
0.031 < T < 0.61 and 0.04 < Ax < 0.43 as the tem- 
perature varies from 4.2 K to 9.2 K. In Fig. ||, T < 0.2 
corresponds to T < 6.7 K and X± < 0.05. At these low 
temperatures, there is a larger variation in T because of 
the sensitive dependence of this parameter to the junc- 
tion's resistance below the gap voltage. 

As Fig. |^ shows, the maximum current supported by 
the breather, I + , is almost equal to Idep- A simple cir- 
cuit model gives some physical intuition. Junctions that 
are rotating have some effective resistance while junc- 
tions that are in the zero-voltage state have zero resis- 
tance. In our breather state, the center junction is ro- 
tating. Therefore, when we apply a current to our array 
the current will tend to flow around the rotating junc- 
tion and through the outside junctions that are in the 
zero-voltage state. When these outside junctions reach 
their critical currents they will begin to rotate and the 
breather will disappear. In the simplest case, when we ig- 
nore any circulating Meissner currents, this model yields 
I+/NI CV = (h + l)/(2h + 1) = 0.8. Since I dep is roughly 
NI CV , the depinning current is the upper bound for the 
applied current that the breather can support. 



The instability mechanism that determines i_ in our 
experiments is more difficult to discern. We offer two 
suggestions that are due to the underdamped charac- 
ter of our system. One possibility is via a retrapping 
mechanism similar to that of a single junction. As 
the middle junction rotates, it reaches a point where 
the current drive is not sufficient to support the ro- 
tation and it destabilizes. This physical picture gives 
I-/NI CV = (2h + 2)(4/7r)r. So that for our parameters, 
I- should be 2.5 times larger than I r , as it is approx- 
imately in Fig. |i| A second possible instability mech- 
anism consist of resonances between the characteristic 
frequencies of the breather and the lattice eigenmodes. 
The breather looses energy as it excites the eigenmodes. 
In our experiments, the breather always looses stability 
at voltages close to V g . For our parameter range, V g is 
larger than the voltages for the lattice eigenmodes, thus 
our data seems to favor a retrapping mechanism. Lastly, 
we add that since Idep and consequently I + decreases 
with r, there also seems to be a critical damping where 
the breather will cease to exist. Experimentally we did 
not find a breather for T > 0.2. 

In summary, we have experimentally detected differ- 
ent breather and multi-site breather states in a super- 
conducting Josephson ladder network. By varying the 
external current and temperature we have studied the do- 
main of existence and the instability mechanisms of these 
localized solutions. In addition we have also found, but 
not discussed here, breathers which are not top-bottom 
voltage symmetric [^2|, in which only the top (bottom) 
horizontal junctions rotate while the bottom (top) junc- 
tions are in the zero- voltage state. These experiments are 
the first observations of discrete breathers and multi-site 
breathers in a condensed matter system. 
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